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^ ■ ABSTRACT 
O 

Q-i! We discuss head-on collisions of neutron stars and disks of dust ("galaxies") 

}L{ '. following the ideas of equilibrium thermodynamics, which compares equilibrium 

states and avoids the description of the djTiamical transition processes between 



them. As an always present damping mechanism, gravitational emission results 
in final equilibrium states after the collision. In this paper we calculate selected 
5^ ! final configurations from initial data of colliding stars and disks by making use 

of conservation laws and solving the Einstein equations. Comparing initial and 
final states, we can decide for which initial parameters two colliding neutron stars 
(non-rotating Fermi gas models) merge into a single neutron star and two rigidly 
rotating disks form again a final (differentially rotating) disk of dust. For the 
neutron star collision we find a maximal energy loss due to outgoing gravita- 
tional radiation of 2.3% of the initial mass while the corresponding efficiency for 
colliding disks has the much larger limit of 23.8%. 



Subject headings: Equation of state — gravitation — gravitational waves — 
galaxies: general — stars: neutron 
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1. Introduction 

Collisions of compact objects are an important source of gravitational radiation. Much 
effort has recently been made to develop numerical methods and codes describing and sim- 
ulating the underlying hydrodynamical and gravitational phenomena. After the pioneering 
work on numerical black hole evolutions by Eppley and Smarr in the 1970's [see e.g. Ep- 
pley (1975) and Smarr et al. (1976)], head-on collisions were re-investigated in the 1990's 
(Anninos et al. 1993, 1995, 1998) with good agreement between numerical and perturbation- 
theoretical results. Long-term-stable evolutions of black hole and neutron star collisions were 
successfully performed in the last two years (Sperhake et al. 2005; Fiske et al. 2005; Zlo- 
chower et al. 2005; Sperhake 2006; Loffler et al. 2006). 

From a mathematical point of view collision processes are typical examples for initial- 
boundary problems. In particular, we will discuss head-on collisions of spheres and disk^, 
see Fig. [H Starting with bodies separated by a large ( "infinite" ) distance we may model the 
initial situation by a quasi-equilibrium configuration of two isolated bodies. Corresponding 
solutions for spheres and (rigidly rotating) disks can be found in the literature, see e.g. Mis- 
ner et al. (2002), Shapiro & Teukolsky (1983) and Neugebauer & Meinel (1993, 1994, 1995). 
The dynamical phase of the collision process is always accompanied by gravitational radia- 
tion. This damping mechanism results again in the formation of an equilibrium configuration 
after the collision. The rigorous mathematical description of the dynamical transition phase 
is difficult and requires extensive numerical investigations. However, interesting information 
about the collision can be obtained by comparing the initial and final (equilibrium) states. 
This thermodynamic idea avoids the analysis of the transition process and reduces the math- 
ematical effort to solving the Einstein equations for the end products, which are stationary 
and axisymmetric in our case. The solution makes use of conservation laws which transfer 
data extracted from the initial configurations (spheres and disks before the collision) to the 
final configurations. 

While the initial configurations are available the calculation of the final states is rather 
difficult. To cope with this problem for head-on colliding stars and disks, we will make use 
of two heuristic principles: 

1) Perfect fiuid configurations at rest are spherically symmetric. Hence, the end product 
of two head-on colliding spheres without angular momentum is again a sphere (a fiuid 
ball or a Schwarzschild black hole). 



^Disk-like matter configurations play an important role in astrophysics, e.g. as models for galaxies, 
accretion disks or intermediate phases in the merger process of two neutron stars. 
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2) Dust configurations are two-dimensional ("extremely flattened") and axisymmetric 
(with non- vanishing angular momentum). Consequently, the dust matter after a head- 
on collision of two disks of dust is again two-dimensional and axisymmetric (a compact 
disk, a disk surrounded by dust rings or a black hole surrounded by dust rings). 

Though plausible, these principles have not been proved rigorously so fail. For proofs under 
special assumptions see Beig & Simon (1992) and Lindblom & Masood-ul-Alam (1994). 

As illustrated in Fig. [1] we will confine ourselves to two problems: 



a) head-on collisions of two identical spheres (stars) merging into a single fluid ball and 

b) head-on collisions of two identical disks of dust (galaxies) merging into a single disk. 



We will be able to formulate necessary conditions for the formation of these balls or 
disks. Obviously, the conditions will restrict the parameters of the initial configuration; 
a violation of the conditions would necessarily lead to other final states such as to black 
holes or central disks surrounded by rings. To express the parameters of the admissible 
initial parameters — the first goal of this paper — we have to solve the Einstein equations 
(numerically but) rigorously and to make use of the conservation laws for baryonic mass 
and angular momentum. There is no obstacle to an extension of the method. One could 
start a systematic investigation of other possible final states after the collision (black holes, 
black holes with rings etc.) making use of symmetries, conservation laws and the heuristic 
principles 1) and 2). An important point of the procedure would be the stability analysis 
of the end products. As for our investigation, there is important evidence from Newtonian 
gravity that rigidly or differentially rotating disks of dust are unstable. Nevertheless we can 



^In this context we refer to an new approach by Masood-ul-Alam (2007). 




Fig. 1. — Model: collisions of spherically symmetric stars or rigidly rotating disks of dust 
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expect that "stabilizing" effects (pressure due to internal kinetic energy) do not falsify our 
other goal — to estimate the maximal contribution of gravitational radiation to the total 
energy loss AM. In general, the total energy loss calculated via the comparison of initial 
and final equilibrium configurations is only an upper limit for the energy loss (efficiency) 
due to gravitational emission. (It includes energy loss due to non-gravitative radiation or 
mass ejection during the collision). We will present a disk collision model which is exclusively 
damped by gravitational radiation. The resulting differentially rotating disk will be compared 
with a rigidly rotating disk of the same baryonic mass and angular momentum formed from 
the same initial disks under the additional influence of dissipative processes in the matter 
(see 13.21) . Thus we can compare the efficiencies of the two forms of dissipation. 



In Sec. [2] we discuss, as introductory examples, the merger of two Schwarzschild stars 
and the collision of two (Fermi gas) neutron stars. Sec. [3] contains the main part of this paper 
which is dedicated to the investigation of disk collisions. These discussions are based on a 
novel solution of the Einstein equa tions describing the final con figuration. Here we continue 



the analysis of a previous paper (IHennig fc Neugebauer 20061 ). in which we discussed the 



collisions of rigidly rotating disks of dust with parallel (or antiparallel) angular momenta 
under the simplifying assumption that the final disk be again a rigidly rotating (or rigidly 
counterrotating) disk of dust. This assumption can only be justified if friction processes 
between the disk rings provide for a constant angular velocity throughout the disk. This 
model seems to be somewhat artificial and unsuited to determining the contribution of 
gravitational radiation to the total energy loss. Interestingly, our present investigation will 
show that the frictional contribution to the total energy loss for colliding rigidly rotating 
disks is comparably small. 



2. Star collisions 

2.1. Introductory example: Schwarzschild stars 

In order to demonstrate the method, we study the collision of two Schwarzschild stars, 
i.e. spherically symmetric perfect fluid stars with a constant mass density, /i = constant. 
Though not very realistic, this model illustrates the main steps of the method. 



The matter of a Schwarzschild star is described by the perfect fluid energy-momentum 
tensor 

= + p)u'u^ + pg'^ (1) 
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with the pressure 



Pir) = \ ^ ^ ^n, (2) 

3^1 - '-frl - y/l - 

where m*, r and tq are the four-velocity, the radial coordinate and the coordinate radius of 
the star, respectively. The interior Schwarzschild metric can be written as 



ds' = + r\d^' + sin^^d^^) -{l^^-^-fri-yi- ^r^j dt^ (3) 

and the exterior Schwarzschild solution is 

= + [M^ + sin" M^p^) - (^1 - dt\ (4) 

Note that we use the normalized units where c = 1 for the speed of light and G = 1 for 
Newton's gravitational constant. 

The gravitational mass M, 

M = ^rl, (5) 

follows from the matching condition at the star's surface and the baryonic mass Mq is given 
by 

Mo= fiu^^/^drdM(p = Anfi / , ^ ^ (6) 

J J Ji-§mr^ 

t=to Y 3 ' 

Now we apply these formulae to the head-on collision of two stars. Restricting ourselves 
to collisions of two identical Schwarzschild stars we assume, that the final star be again a 
Schwarzschild star and have the same mass density (e.g. nuclear matter density), 

fi = fi, (7) 

where from now on tildes denote quantities after the collision. 

The conservation of baryonic mass during the collision process 

Mo = 2Mo, (8) 

allows one to calculate the parameters of the final star as a function of the initial parameters. 
With (IS]), ([6]) and (I7l) the conservation equation ([H]) can be written as 



, /2Mfo\ /2Mfo 2Mfl / /2M /2M / 2M , , , 
arcsin \ ^ - J ^Wl 1 = 2 arcsin . / J Jl , 9 
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i.e. the radius ratio ro/ro is a function of the initial mass-radius ratio 2M/ro. Hence, we 
may express the efficiency rj of conversion of mass into gravitational radiation, 

,.i_4.i4ff.V. (10) 



2M 2 \ro 

and the mass-radius ratio of the final star, 

H^ = H^iff^y. (11) 

^"0 ro \roJ 

in terms of 2M/rQ. 

The resulting parameter relations are plotted in Fig. [2l For Schwarzschild stars the 
coordinate radius is restricted by the Buchdahl condition, i.e. 

ro > ^ X 2M, ro > ^ X 2M. (12) 

o o 

As a consequence, the first plot shows, that "relativistic" initial stars with 2M/ro > 0.6482 . . . 
can never merge into a new Schwarzschild star with the same mass density /x. The "physical" 
parts of the parameter relations are shown as solid curves while the forbidden parts are 
dashed. According to the third plot the efficiency rj cannot exceed a maximal value of 
Vmax ~ 19.7%. 



2.2. Neutron stars: Completely degenerate ideal Fermi gas 

In order to extend the discussion of the previous section to a more realistic star model, 
we replace the equation of state fi = constant by the equation for a completely degenerate 
ideal Fermi gas of neutrons. 

The (interior) line element of a spherically symmetric star can be written as 

ds^ = e^^^^^dr^ + r\di}^ + sin^ ^ dcp^) - e^'^'^df (13) 
and the matter is again described by the perfect fiuid energy-momentum tensor 

= {fi + p)uV + pg'^ . (14) 
With the definition of a new metric function m{r) by 
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2M/ro 

Fig. 2. — Parameter relations for colliding Schwarzschild stars: The final mass-radius ratio 
2M/fo, the radius ratio fo/f and the efficiency r) are plotted as functions of the initial 
mass-radius ratio 2M/rQ. Dashed parts of the curves mark regions inaccessible due to the 
Buchdahl inequality 2M/fo < 8/9. 
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the field equations can be written in the TOV form, see e.g. Shapiro & Teukolsky (1983), 

— = 47rrV, "^(0) = (16) 
dr 

dp m f p\ f 47rr^p\ / 2m\ ^ , , , _, 



dr 



where pc is the pressure in the center of the star. 

We will solve these equations for the completely degenerate ideal fermi gas of neutrons 
with the equation of state [see e.g. Shapiro & Teukolsky (1983)] 

p = cif{x), p = c2X^ iJ, = p + cig{x), (19) 

where 

f{x) = x{2x^ - 3) Vl + x2 + 3 ln(x + Vl + x^), (20) 



g{x) = 8x^(^1 + a;2 - 1) - /(x), (21) 

with the neutron mass mn = 1.6749286 x 10^'^ kg and Planck's constant h = 6.626076 x 
10~^^ Js. By solving the TOV equations (fT6l) and (ITTI) with the equation of state (fT9|) for 
a sequence of values of the central density one can calculate the corresponding radii of the 
stars as the first zero tq of p{r), their gravitational mass from M = m(ro), and their baryonic 
mass as 

ro 

Mo = 4n [ _P^t^. (23) 

The resulting mass-radius relations are shown in the first plot of Fig. [31 
Again the baryonic mass is an invariant of the collision, i.e. 

Mo = 2Mo (24) 

for the collision of two identical initial stars. This equation has to be analysed together 
with the mass-radius relations. (Thereby we take into account only stars in the monotonic 
decreasing part of the mass-radius relation Mo(ro).) The resulting parameter relations are 
shown in the remaining plots of Fig. [3l For the maximum of the efficiency one finds r/max ~ 
2.3%, i.e. a comparably small value in view of the maximal efficiency ?7max ~ 19.7% for 
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2M/ro 2M/ro 



Fig. 3. — Parameter relations for the collisions of neutron stars made up of degenerate 
neutrons [cf. f|T9l) ]. First plot: mass-radius relations for the baryonic mass Mq and the 
gravitational mass M. Second plot: initial radius tq and final radius as functions of 
the mass-radius ratio 2M/rQ. Third plot: change of the coordinate radius. Fourth plot: 
efficiency rj compared to the efficiency of the collision of Schwarzschild stars (dashed curve). 
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collisions of Scliwarzschild stars. The reason is the relatively small maximal mass of Mmax ~ 
0.7 Mq permitted by the equation of state f|T9|) that excludes highly relativistic values for 
the mass-radius ratio 2M/rQ. However, compared to Schwarzschild stars with the same 
parameter 2M/ro the collisions of neutron stars are more efficient, cf. the last plot in Fig. [31 

Another difference is the change of the coordinate radii. While two Schwarzschild stars 
merge into a Schwarzschild star with a coordinate radius bigger than the initial radius, 
'f^o/fo > 1 (cf- Fig. [2]), the resulting neutron star is smaller than the initial neutron stars, 
fo/ro < 1 (cf. Fig. ED. 

3. Disk collisions 

Collisions of disks of dust require more effort. In particular, the discussion of the final 
equilibrium state is based on a solution of a free boundary value problem to the Einstein equa- 
tions. At the first glance, the conservation laws for baryonic mass and angular momentum 
are not sufficient to formulate a complete set of boundary conditions for the configuration 
after the collision. However, excluding non-gravitational dissipation, we may replace the 
global conservation laws, as used in Sec. [21 by local ones. Due to the geodesic motion of dust 
particles, the baryonic mass and the angular momentum of each of the rings forming the disk 
are conserved separately, see Fig. [3 Using such local conservation laws we will be able to 
solve (numerically) the boundary value problem for the final state after the head-on collision 
of two aligned rigidly rotating disks of dust with parallel angular momenta, cf. Fig. [H 

3.1. Initial disks: Rigidly rotating disks of dust 

The free boundary value problem for the relativistic rigidly rotating disk of dust (RR 
disk) was discussed by Bardeen and Wagoner (1969, 1971) using approximation methods 
and analytically solved in terms of ultraelliptic theta functions by Neugebauer and Meinel 
(1993, 1994, 1995) using the Inverse Scattering Method. The line element of the stationary 
(Killing vector: ^*) and axisymmetric (Killing vector: rj^) space-time may be written in the 
Weyl-Lewis-Papapetrou standard form 

ds2 = e-2f^[e2'=(dp' + dC') +pW] -e'^(dt + ady^)2, C = ^1 = ^ip, (25) 

where the metric potentials U = U{p,(), k = k{p,() and a = a{p,() are given in terms of 
ultraelliptic theta functions. 
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The matter of the disk of dust is described by the energy-momentum tensor 

= e{p)6{C)u'u\ (26) 

where e{p)S{() is the mass density with 6{() as Dirac's (5-distribution. Due to the symmetries, 
the four-velocity of the dust particles is a linear combination of the two killing vectors, 

= e-^°(f + nor]'), u'ui = -1, (27) 

whence 

(e + nor?')te + ^^or?.) = -e'''°, (28) 

where VLq is the angular velocity of the particles forming the disk and Vq is a redshift param- 
eter. Rigid rotation means VLq = constant in the disk. Since dust particles move geodesically 
this assumption implies Vq = constant in the disk. Hence, the boundary condition fl28l) and 
as a consequence the RR disk solution contains two constant parameters. Alternatively to 
Qq and Vq, one may choose the coordinate radius po of the disk and a centrifugal parameter 
/i = 2Qlple~'^^° [/i — >• turns out to be the Newtonian limit and /i — >• 4.62966 . . . the ultra- 
relativistic limit, where the disk approaches the extreme Kerr black hole, cf. Neugebauer & 
Meinel (1994) and Neugebauer et al. (1996) for these and further properties]. 



3.2. Final disk: Differentially rotating disk of dust 



In a previous paper (IHennig fc Neugebauer 20061 ) we discussed head-on collisions of two 



(identical) rigidly rotating disks of dust merging into one rigidly rotating disk of dust. The 
model excluded mass ejection and made use of the conservation of baryonic mass and angular 
momentum (axisymmetry). From a thermodynamic point of view rigid rotation of the final 
disk means thermodynamic equilibrium, which is a result of dissipative processes during 
the dynamical phase. Hence, the energy difference between the initial state (two separated 
disks) and the final state (one rigidly rotating disk) is influenced by irreversible processes 
in the matter and outgoing electromagnetic radiation as well as by emission of gravitational 
waves. The intention of this paper is to compare the contribution of these two effects by 
calculating the end product of a purely gravitational collision process which we expect to 
be a differentially rotating disk of dust. Note that our thermodynamic analysis enables us 
to formulate necessary conditions for the parameters of the initial disks (/x restricted) to 
permit the formation of a final disk. To obtain sufficient conditions one would have to solve 
the Einstein equations for the time-dependent transition phase, which is outside the scope 
of this paper. 
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In the next subsection we will give a brief summary of the previous paper. After that, 
we will see that the local conservation of baryonic mass and angular momentum is sufficient 
to calculate the final differentially rotating disk (nur aerically) . Differentially rotating disk s 
with arbitrary rotation law have already been studied (lAnsorg fc Meinel 200d : lAnsorg 200ll ). 
The point made here is that we are able to formulate a physically motivated rotation law as 
a result of a collision process. 



3.2.1. Formation of rigidly rotating disks 

For the formation of an RR disk from two colliding RR disks under the influence of 
friction processes the conservation equations for baryonic mass and angular momentum, 

Mo = 2Mo, J = 2J, (29) 

are sufficient to calculate the parameters of the final disks as functions of the initial pa- 
rameters. These equations and explicit formulae connecting the gravitational mass M, the 
baryonic mass Mq and the angular momentum J of the RR disk allowed us to calculate 
the efficiency 

^rvrv ^ i_ M/2M as a function of the initial centrifugal parameter /z, cf. 
Fig. m It should be emphasized once again, that this efficiency measures the total energy 
loss including friction. Therefore rj is only an upper limit f or the energy of the gravitati onal 



emission. We obtained a maximal value of ^/max 

^ 23.8% flHennig fc Neugebauer 20061 ). 



Furthermore, it turned out that the formation of RR disks from two colliding RR disks 
is only possible for a rather restricted interval < /i < 1.954. . . of the initial centrifugal 
parameter /i. If /x exceeds this limit, the collision must lead to other final states, e.g. black 
holes or black holes surrounded by matter rings. 



3.2.2. Local conservation equations 

We now turn to the main goal of this paper and analyse the formation of a disk of 
dust under the influence of gravitational forces as the only form of interaction. Comparing 
the resulting differentially rotating disk of dust with the rigidly rotating disk of the same 
baryonic mass Mq and angular momentum J formed from the same initial disks we may 
separate gravitational damping due to the emission of gravitational waves from frictional 
processes in the matter. 

We may interpret a disk of dust as a superposition of infinitesimally thin dust rings. 
Considering the geodesic motion of a single mass element, one can show that for correspond- 
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Fig. 4. — The efficiency jyj^j^ for the formation of an RR disk from two initial RR disks as 
a function of the centrifugal parameter of the initial disks, rj ^ is an upper limit for the 
energy loss due to gravitational radiation. 




Fig. 5. — Illustration of the local conservation equations. Two corresponding rings of the 
initial disks merge into a ring in the final disk. The baryonic mass dMo and the angular 
momentum d J of these rings are conserved. 
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ing rings in the two initial disks (see Fig. [5]) the baryonic mass and the angular momentum 
are conserved, 

dMo = 2dM, dJ = 2dJ, (30) 

i.e. the baryonic masses dMo and the angular momenta dJ of the rings with radius p, 
taking up the interval [p, p + dp], in each of the two initial disks sum up to dMo = 2dMo 
and dJ = 2dJ of the corresponding ring in the final disk (with radius p, taking up the 
interval [p, p + dp]), cf. Fig. O It should be emphasized that the local conservation laws 
f l30|l would be violated by dissipative processes in the matter or, mathematically speaking, 
by dissipative terms in the total energy-momentum tensor as the source of the Einstein 
equations during the collision phase. Having reached a final equilibrium configuration (e.g. 
a rigidly rotating disk of dust) the system "forgets" the dissipative terms and behaves like 
cold dust with an energy momentum tensor of the form (!26|) . During the interaction phase, 
angular momentum will be transported within the disk by viscous forces and only the total 
angular momentum (axisymmetry!) and the total baryonic mass are conserved (!29|) . The 
ring-wise conservation of baryonic mass and angular momentum fl3UI) is characteristic for 
purely gravitational damping processes. They arise from collision processes governed by an 
energy-momentum tensor of dust without dissipative terms. In this case the geodesic motion 
of the volume elements implies the conservation of baryonic mass and angular momentum 
in each volume element and therefore implies fl30|) . 

Eq. ( 130|) provides us with a subset of the boundary conditions to be discussed in the 
next subsection. It will turn out that these conditions, together with conditions resulting 
from the field equations, determine a unique solution of the Einstein equations describing a 
final disk with differential rotation (DR disk) as the end product of the collision process. 



3.2.3. Boundary value problem for the final DR disk 

The line element f l2B]) . which may also be used to describe any axisymmetric and sta- 
tionary differentially rotating disk, can be reformulated to give 

d~s^ = e2^(dp2 + dC^) + p2e-2^(d<^ - Cjdif - e'^dP, (31) 

where the usage of the functions k, v and u; [instead of \J , k and a as in fl25l) ] avoids numerical 
issues with ergospheres (where e^^ < 0). According to ( l26l) the energy-momentum tensor is 

= e{p)6{C)u'u^ (32) 

and the four-velocity is again [cf. (1271) ] a linear combination of the killing vectors, 

u' = e-^{e + Qv'), (33) 
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where V = V{p) and Q = ^l{p) are functions of p [constancy of V and Q defines rigid 
rotation, cf. (|271)]FI 

The vacuum field equations for u and u are [cf. Bardeen (1973)] 

Ai z> = je-^'iCoJ + Cop, AscD = 4(9 pCo^^ + u^^u^^), (34) 

with 



71 

A„ := dl + dl + -rd^. (35) 
The matter appears only in the boundary conditions along the disk {( = 0, p < po), 

^,clc=o+ = -8^^r^' (37) 
p{Q - = (1 + v'^y^^v-p +~p^{Sl- Cj)Cj^p, (38) 

where 

v:=pe-^''{VL~Cj), a := ee^~\ (39) 
Thus we have to deal with a boundary value problem for Einstein's vacuum equations. 

As already mentioned, the local conservation equations (l30l) of the previous subsection 
lead to additional boundary conditions along the disk. From dMo = 2dMo with dMo = 
27r(je~^°pdp and dMo = 27Tae~^ pdp we obtain 

pe^°dp , 
a = 2a^-^. 40 
pe^dp ^ ' 

Likewise, dJ = 2dJ with dJ = 27i(7e~^°u^r]ipdp and dJ = 27rcre~^£t*f/jpdp, w^r/j = pfe"^" and 
u^Vi = pve~^ leads to 

pye"^ = pwe"^". (41) 
The function V^(p) can be calculated from -u^-Uj = —1, 

e^y = (1 - {)2)e2P_ ^42) 

The remaining boundary conditions describe the behaviour at infinity, where the metric 
approaches the fiat Minkowski metric, 

k = i> = u = 0, (43) 



■^AU quantities of the final DR disk are tilded. 
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and in the plane C = outside the disk (p > po), where fl36l) and fl37|) lead to vanishing 
normal derivatives, 

^.dc-o.=0' ^.dc-o.=0- (44) 
In addition we have to ensure regularity along the axis of symmetry p = 0. 

Eqs. (I3ll) . (I36l) - fl38l) and fj40|) . fj4T|) form a complete set of equations to determine the 
unknown functions uniquely: There are two two-dimensional functions, I'ipX) '^(p, C); 
which have to satisfy the two elliptic partial differential equations flM|) with the boundary 
conditions and fl57|) . and three additional one-dimensional functions in the disk, f2(p), 
cr(p), p(p), which have to obey the three boundary conditions (138|) . (HOj) and (14T]) . (The 
metric function k can be calculated by a line integral afterwards, but is not needed for the 
computation of the efficiency t] in our collision scenario.) 



3.2.4- Numerical method 



In order to prepare numerical investigations we will map the region < p < oo, < 
C < C)0 to a unit square thus reaching a compactification of infinity, cf. Fig. [61 (Due to the 
refiection symmetry with respect to the plane ^ = we can restrict ourselves to the region 
^ > 0.) To do this we introduce in a first step elliptical coordinates 

p=V{'^ + e)i'^-v'), C = ^V, ee[0,oo], r/G [0,1] (45) 

(without loss of generality, we may choose units where po = 1). In a second step we stretch 
the coordinates by the transformation 

e = cot (^s) , r] = VT^, s G [0, 1], t G [0, 1]. (46) 



disk £ p 




oo 



Fig. 6. — The coordinate transformation fl45l) . fl46|) maps the part C ^ of the p-^-plane to a 
unit square in the s-t-plane. £ denotes the equatorial plane outside the disk, ( = 0, p > po. 
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The coordinates s and t form a unit square with the following boundaries, 
s = : oo 

s = 1 : disk, p < 1, C = 

t = : axis of symmetry, p = 

t = 1 : disk plane £ outside the matter, p > 1, C = 0. 

The unknown functions in the boundary value problem are analytic functions in this square 
(as is known for the case of Maclaurin disks or the RR disks). Hence, it is convenient to use 
spectral methods for the numerical solution of the boundary value problem. We expand the 
unknown potentials in terms of Chebyshev polynomials Tj to a predetermined order in the 
form 

f{s,t) = J2cjkTj{2s-l)Tk{2t-l) or /(t) = ^ CfcTfc(2t - 1) (boundary) (48) 

j,k k 

and formulate the Einstein equations at the extrema of the Chebyshev polynomials. This 
leads to an algebraic system of equations for the Chebyshev coefficients (or, alternatively, 
for the values of the potentials at these points) that can be solved with the Newton-Raphson 
method. The iteration starts with an initial "guessed" solution (for example the Newtonian 
approximation, see Sec. 13.2.61 below). 

The calculations show a decreasing accuracy of the numerical solution for increasingly 
large values of the initial parameter p. The reason are large gradients of the metric potentials 
for strong relativistic DR disk which make the Chebyshev approximation more costly. To 
reach a better convergence we perform an additional coordinate transformation 

^ = (49) 

smh(oj 

introducing a new coordinate s, where 5 is a suitably chosen parameter. As shown in 
Ansorg & Petroff (2005), this transformation smooths the gradients of the metric functions. 
The convergence is illustrated in Fig. [71 



3.2.5. Results 

Using this numerical algorithm, we are able to solve the boundary value problem for the 
final DR disk. In particular, we could calculate, for each value of the initial parameter p, 
all metric coefficients of this final disk. However, we will restrict ourselves to the discussion 
of the relations between the initial and final parameters and the efficiency of the collision 
process. Especially, we will compare the final DR disk with an RR disk having the same 
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Fig. 7. — Convergence properties of the numerical code for the example rj . The values of 

the efficiency rj for different orders Ug = rit = n of the Chebyshev expansion are related 
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to the order n = 32. The plot shows 
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Fig. 8. — Two models for disk collisions: 

(A) Under the influence of a small amount of friction, RR disks merge again into an RR 
disk. This scenario was discussed in Hennig & Neugebauer (2006), see Sec. 13.2.11 

(B) In the absence of friction, the same RR disks merge into a DR disk. Allowing for friction 
afterwards, the system would again arrive at the RR disk of scenario (A) after a sufficiently 
long time. 
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baryonic mass and angular momentum. The point made here is that such a rigidly rotating 
disk represents the state of "thermodynamic equilibrium" for disks of dust as the end point 
of their thermodynamic evolution. As sketched in Fig. [HI there are at least two possibilities 
for the formation of this final RR disk: The direct process (A) including friction from the 
beginning or the equivalent thermodynamic process (B) where, in a first step, a differentially 
rotating disk is formed (by gravitational damping alone, no friction) and, in a second step, 
the angular velocity becomes constant (due to friction). Note that baryonic mass and angular 
momentum are conserved in both processes. By comparing (A) and (B) we may extract the 
contribution of friction in scenario (A). 

In the following discussion, tilded quantities, as before, belong to the final DR disk of 
scenario (B) in Fig. [HI a superscript "RR" denotes quantities of the final RR disk in scenario 
(A) and the centrifugal parameter /i without any additions characterizes the initial RR disks. 

The rotation curve of the final DR disk, i.e. its (normalized) angular velocity Qpo 
as a function of the (normalized) radius p/po is shown in the first two plots of Fig. [9] 
For small parameters p (post-Newtonian regime) the function ClpQ is almost constant (first 
plot). Interestingly, strongly relativistic disks (/i > 1.5) show the same property (second 

plot). Moreover, Qpo tends to zero in the ultrarelativistic limit in analogy to the relation 

R R R R R R 

j^rtrt^rtrt _^ g ^j-^jch holds for RR disks in the ultrarelativistic limit 4.62966 . . . 

The "centrifugal parameter" p = 2Cl'^pQe~^ = p{p) is shown in the third plot of Fig. [Hi 

Like the angular velocity, p is almost constant for small p. For strongly relativistic DR disks, 

• • R R 

p in the center of the disk exceeds the limit Pm^xx ~ 

4.62966 . . . of RR disks. 

The fourth plot of Fig. [9] shows the quantity QM as a function of p/po- For strongly 
relativistic DR disks ClM becomes constant and approaches the limit 0.5. On the other 
hand, this is a characteristic value for extreme Kerr black holes where ^H"^BH ~ (^H' 
angular velocity of the horizon, Mgjj: black hole mass). Indee d, one can show that there 



is a phase transition betwe en RR disks and Kerr black holes (iBardeen fc Wagoner 1971 



Neugebauer fc Meinel 19941 ). This inspires the conjecture that the DR disk exhibits the same 
phase transition. There is no obstacle for a (numerical) proof of this assumption in principle. 
To extend our present code to study the parametric collapse of the DR disk including the 
formation of a horizon we would have to follow the ideas of Bardeen and Wagoner (1971) 
who analysed this problem for the RR disks. However, such investigations are outside the 
scope of this paper. 

The fifth plot of Fig. M shows the redshift z for a photon emitted from the disk center 
as a function of the initial centrifugal parameter p. For increasing values of p (relativistic 
DR disks) z grows rapidly. 
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Fig. 9. — Parameter relations for the collision of RR disks. We performed numerical calcu- 
lations for values of the initial centrifugal parameter /i in the invervall [0, 1.9]. The dotted 
part of the curve in the last plot is an extrapolation for larger fi. This extrapolation and 
the rapidly growing redshift in the second last plot indicate, that the initial parameter /i in 
scenario (B) is limited (approximately, or perhaps even exactly) to the same interval as in 
scenario (A), < ;U < 1.954 . . ., see Sec. 13.2.11 
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An important result is the efficiency 77 of the formation of DR disks which measures 
the amount of energy converted into gravitational radiation. The difference 1]^ = j^R-R-— j^DR 

ID ID 

as shown in the last plot of Fig. [9] compares this value with the efficiency of the RR 

disk forming process as sketched in scenario (A) of Fig. [HI Thereby, t]^ is the part of energy 

F R R F)R 

lost due to friction during the formation of a final RR disk. We find r] = rj — fj^^ < 
1.5 X 10~^, i.e. the contribution of friction is extremely small, t]^ <^ ^RR^ such that the 
gravitational radiation dominates the collision process (A). 



3.2.6. Analytical treatment of the Newtonian limit 

Our numerical investigations have shown that the angular velocity of the final DR disk 
becomes closer and closer to a constant over the whole range of p/po as the centrifugal 
parameter /i tends to zero (cf. the ffist plot of Fig. This leads one to suspect that a 
final disk with a strictly constant angular velocity will solve the boundary value problem 
as discussed in Sec. 13.2.31 in Newtonian theory. Interestingly, we can treat this problem 
analytically. This will now be demonstrated. Strictly speaking, there is no gravitational 
radiation in Newton's theory. However, this Newtonian boundary value problem can be 
seen as the limit of a sequence of relativistic collisions with decreasing p, all reaching a final 
equilibrium state due to gravitational emission. Moreover, the Newtonian solution can be 
used as a starting point for the iterative calculation of the final relativistic DR disk. 

Since the Newtonian limit of the RR disk is the Maclaurin disk we have to study 
the collision of two identical Maclaurin disks using the local conservation laws (I5U]) . The 
Newtonian potential U of the final disk is a solution of the Poisson equation 

A t/ = 47ra(5(C) (50) 

with the boundary condition 



= 27r5-, (51) 

C=o+ 



where a = a{p) is the surface mass density of the final disk. With dM = 27rcrdp and 
dJ = Qp^dM, Eq. fl5Ul) leads to the additional boundary conditions 

a(p)=2a(p)?^, Q{p)p' = nop'. (52) 
pdp 

The initial surface mass density of the Maclaurin disk is 

^(p) = 7?^Jl-4 (53) 
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and the initial constant angular velocity is related to the initial mass by 
Using these relations, together with the Euler equation 



Up 



- =^\~P)P, (55) 
C=o 



we find that a (rigidly rotating) Maclaurin disk with the parameters 



n = 4^0, Po = ^Po (56) 



indeed solves the boundary value problem (l50l) - (!52l) . 



4. Discussion 

In this paper we have performed the analysis of collision processes in the spirit of equi- 
librium thermodynamics. Avoiding the solution of the full dynamical problem, we compared 
initial and final equilibrium configurations to obtain a "rough" picture of these processes. 
In this way we were able to calculate the energy loss by the emission of gravitational waves 
and to find conditions ( "parameter relations" ) for the formation of final stars and disks. 

The application of this method to collisions of perfect fiuid stars and collisions of rigidly 
rotating disks of dust leads to restrictions of the initial parameters. It turned out that the 
formation of final stars/disks from stars/disks is only possible for a subset of the parameter 
space of the initial objects. Otherwise, the collision of spheres and disks would lead to other 
final states, e.g. to black holes. 

Our main result is the numerical solution of the Einstein equations for the differentially 
rotating (DR) disk formed by the collision of two identical rigidly rotating (RR) disks with 
parallel angular momenta. We calculated the characteristic quantities of the final DR disk, 
as for example the rotation curve fl{p) as it depends on the centrifugal parameter fi of the 
initial RR disks. It turned out, that the angular velocity Q is almost constant (as shown in 
Sec. 13.2.61 it is strictly constant in the Newtonian limit). Therefore, the simplified model of 
the formation of an RR disk from the collision of two RR disks as presented in Hennig & 
Neugebauer (2006), which has to allow frictional processes to reach constant angular velocity, 
turns out to be a good approximation to our present purely gravitational (frictionless) model 
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colliding objects 


'7max 


Schwarzschild BHs 


29.3% 


RR disks 


23.8% 


Schwarzschild stars 


19.7% 


Neutron stars 


2.3% 



Table 1: Upper limits for the efficiency t] of different collision processes including 
Hawkins's and Ellis' limit for the collision of two spherically symmetric black holes 



f lHawking fc Ellis 19731 ). According to the last plot of Fig. M [r/^(1.954. . .) = 0] the two 



efficiencies rj^^ and r]^^ coincide with a maximum value ?7max — Virmx ~ 23.8%. 



For each of the studied collision scenarios, we calculated an upper limit for the energy of 
the emitted gravitational waves. A summary of the maximal efficiencies is given in table [TJ 
The value //max ~ 2.3% for the collision of Neutron stars is relatively small compared to the 
other examples. The reason is the restricted equation of state (completely degenerate ideal 
Fermi gas) that does not allow for strongly relativistic stars. 



We would like to thank David Petroff for many valuable discussions. This work was 
supported by the Deutsche Forschungsgemeinschaft (DFG) through the SFB/TR7 "Gravi- 
tationswellenastronomie" . 



A. Potentials of the rigidly rotating disk of dust 

For the numerical calculation of the DR disk that is formed by the collision of two RR 
disks we need some formulae for quantities of the RR disk of dust. 



The coefficient Vq in the four-velocity (l27j) as a function of the p arameter /i can be 
calculated from a very rapidly converging series, cf. (IKleinwachter 19951 ). 

coth ^ = -1 + 0.0294938052100425142/i + 5.4681333461446 ■ 10~V^ 
2 /i 

-1.07467432587 ■ lO^V^ + 2.1127368 ■ 10~^V (Al) 
-4.154- 10-^V + C(/^)- 

The disk values {( = 0, p < po) of the metric functions U and a and the mass density a are 
given by the equations 

= e'''°^f'^-^, (A2) 
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with 



^^0 



'1 + Q,oa)e 



2U 



eVb(/^)gVb(A) 



a 



Vo 



^0 b'oifi) 



/i e 
2V0 



Any, 



(A3) 
(A4) 

(A5) 



cf. (INeugebauer fc Meinel 1994J ). The notation Vo(/i), &o(/^) indicates that the argument fi 
in the parameter functions Vo(/i) and 6o(/i) has to be replaced by /i = (1 — p'^/pDfi. &o(/i) 
means d6o(/i)/d/i. 
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